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Abstract 

We propose a novel way to formulate Chiral Perturbation Theory in 
a nuclear background, characterized by a static, non-uniform distribution 
of the baryon number that describes the finite nucleus. In the limiting 
case of a uniform distribution, the theory reduces to the well-known zero- 
temperature in-medium ChPT. The proposed approach is used to calculate 
the self-energy of the charged pion in the background of the heavy nucleus 
at 0{p^) in the chiral expansion, and to derive the leading terms of the 
pion-nucleus optical potential. 
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1 Introduction 



The topic of low-energy pion-nucleus interactions has a long history and a well- 
developed phenomenology IH El Ej- Energy spectra of pionic atoms have tradi- 
tionally played an important role in promoting the understanding of this field by 
providing a large and systematic data base which sets tight constraints on the 
parameterization of the pion-nuclear potential (see Ref. |1] for a state-of-the-art 
update). 

New developments, both experimentally and theoretically, have recently revi- 
talized this field. Accurate measurements of deeply bound (Is) states of pionic 
atoms formed with Pb and Sn isotopes [S] have triggered renewed interest in the 
underlying mechanisms governing S'-wave pion-nucleus interactions. The quest 
for "fingerprints of chiral restoration" in this context has been raised |H1IZI, re- 
ferring to a possible in-medium reduction of the pion decay constant which is an 
order parameter of the spontaneously broken chiral symmetry in QCD [Hj. Re- 
cent theoretical investigations OEIEIIH based on Chiral Perturbation Theory 
(ChPT) have, for one part, focused on the calculation of the in- medium shift of 
the pion mass, to be identified with the value of the threshold S'-wave vr-nucleus 
optical potential in the center of the nucleus. Alternatively, the energy- dependent 
pion self-energy (polarization function) has been extrapolated, using ChPT input 
and the local density approximation, to calculate directly the pionic Is and 2p 
level shifts and widths for Pb and Sn isotopes P^ . 

What all these recent investigations have in common, is their emphasis on 
a framework based on chiral dynamics and on the spontaneously broken chiral 
symmetry of low-energy QCD. However, all these analyses rely, in one way or 
another, on assumptions which need further systematic scrutiny and clarification: 

i) The in-medium mass of the pion is a well-defined quantity within the 
in-medium ChPT. The procedure of extracting the empirical value for 
this quantity, however, involves detailed theoretical analysis based on phe- 
nomenological potentials, introducing model dependence which is hard to 
control. 

ii) The wave function of a deeply bound pionic state is concentrated near the 
surface of the nucleus, where the variation of the proton and fermion density 
distributions is maximal jTH]. The commonly used expansion in gradients 
of the local nuclear density must therefore be carefully examined. 

iii) Pion-nuclear observables should not depend on the parameterization of the 
interpolating pion field in the chiral Lagrangian. Yet, questions concerning 
off-shell ambiguity have appeared repeatedly in recent works. In Ref. ^2] 
it is demonstrated that the off-shell ambiguity disappears in the mass (as 
it should of course), if one considers the systematic chiral expansion of the 
in-medium mass shift. It remains to prove the equivalent statement for the 
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bound-state mass spectrum, e.g. as obtained in the approach of Ref. [IS]. 
Last but not least, the complete expression of the pion self-energy opera- 
tor at O(p^) and 0{p^), including all strong and electromagnetic isospin- 
breaking effects at this orders, cannot be found in the literature so far. 

It is therefore evident that further progress in the description of pion-nuclear 
bound systems requires the formulation of the systematic ChPT framework in 
the presence of a finite nucleus. This is what our present paper is focused on. 
The interest in this problem is motivated by at least three reasons: 

i) It has been argued (e.g. jlSj) that chiral dynamics (the specific energy 
dependence of low-energy pion-nucleus interactions as imposed by chiral 
symmetry, in combination with the approximate vanishing of the isospin- 
even ttN amplitude) is a key ingredient in understanding the empirically 
observed "missing repulsion" in the S'-wave pion-nucleus optical potential. 
This statement has to be verified on a more systematic basis. 

ii) The in-medium mass shift of the charged pion is the only quantity known 
to us, which is related to the two-point Green function and where the O(p^) 
ttN electromagnetic low-energy constant (LEG) /i enters at lowest order 
independently from other unknown LEGs. On the other hand, this constant 
is very important in the theoretical analysis ^Tj of the experimental results 
obtained by the Pionic Hydrogen GoUaboration (PSI) JH] (in particular, it 
provides the bulk of uncertainty in the relation of the ground-state energy- 
level shift to the S'-wave vrA^ scattering lengths [IZ|). One can expect that 
reliable experimental information on the in- medium mass shift can, at least, 
test the order-of- magnitude estimates for /i, on which the current numerical 
analysis of pionic hydrogen is based |T71 IT^ 1201 • 

iii) The consistent formulation of GhPT in the background of a finite nucleus 
is a challenging task of its own. It opens perspectives for many applications 
of GhPT in nuclear physics - especially for the problems in which the finite 
boundaries of the system are crucial. 

The problem that we are addressing can be formulated as follows. At low 
energy the interactions between pions, photons and nucleons are described by 
the conventional Lagrangian of GhPT (we concentrate here exclusively on the 
2- flavor case). Let Q be the bound state of a nucleus containing A nucleons. We 
aim to describe the processes Q + X ^ Q' + Y where X, Y represent any number 
of pions and photons (in this paper we consider the case where both X and Y 
stand for 1-pion states). In principle one can generalize the approach to cases 
such as pion inelastic scattering or pion absorption, where unbound nucleons are 
also present. This is, however, beyond the scope of the present paper. 
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A rigorous field theoretical treatment would require solving the nuclear bound- 
state (e.g. Bethe-Salpeter) equation and then calculating transition matrix el- 
ements with in- and out-going pion fields. This task is simplified considerably 
for heavy nuclei (in the large- A limit). In particular, in this limit the recoil of 
the nucleus as a whole is a 1/A correction and can be neglected, and one may 
formally set Q' = Q. ChPT in a nuclear background which will be referred to as 
"Chiral Perturbation Theory for Heavy Nuclei" (ChPT® , where @ denotes the 
heavy nucleus with the mass number A) hereafter, is then to be understood as an 
approximate theory in which the transitions Q + X — > Q + Y are treated in terms 
of simpler Green functions "in the presence of the nucleus", Gn{X Y). The 
process X ^ Y is described by chiral dynamics, and the presence of the nucleus Q 
is parametrized in terms of a few phenomenologically determined functions (e.g. 
proton and neutron distributions) which are directly related to corresponding 
observables. 

In-medium ChPT for a homogeneous system starts from a nuclear Fermi 
gas and systematically introduces interactions mediated by real or virtual pi- 
ous. ChPT® is a generalization for finite systems. It reduces to the standard 
in-medium ChPT in the limit of a uniform distribution of baryon number. Be- 
yond that, ChPT® provid es a systematic framework to study pion-nuclear bound 
states, for which the finite volume and the surface of the nucleus are important 
ingredients. 

ChPT® imposes chiral counting rules not only on the hierarchy of pion- 
nucleon interactions, but also on the relevant nuclear matrix elements. The 
nuclear structure information required is thus limited to expectation values of 
only those operators whose chiral dimension is compatible with the given order 
in the chiral counting at which the calculations are carried out. We demonstrate 
explicitly how this works in the calculation of the pion self-energy at 0{p^) in 
ChPT. To that chiral order, the complete set of leading terms in pion-nucleus 
optical potential (those linear in the proton and neutron densities, Pp(r) and 
p„(r)) are generated. It is of course well known from pion-nuclear phenomenol- 
ogy [111210]; that rescattering and absorption terms of 0{p^) and higher orders 
are quantitatively important. In the present paper we focus on the systematic 
expansion to 0{p^). Further developments concerning higher orders are planned 
for subsequent publications. 

The layout of the paper is as follows. In section |2] we consider in detail 
the construction of the framework: main assumptions, role of chiral symmetry 
and the connection to the conventional in-medium ChPT. In section El using 
the formulated framework, we obtain the complete expression for the pion self- 
energy operator at 0{p^) in ChPT in the background of the finite nucleus. Using 
the latter expressions the pion-nucleus optical potential at 0{p^) is derived in 
sectional In section we compare with the existing approaches in the literature. 
Conclusions are drawn in sectional 



4 



2 The framework 



2.1 iS-matrix as functional of the free fields 

In order to ensure that the Green functions satisfy constraints imposed by chiral 
symmetry at every step of the calculation, and do not depend on the param- 
eterization of the pion field, we equip the effective Lagrangian with external 
c-number sources and consider the generating function in the presence of these 
sources. Thus, our approach is closely related to the approach of Refs. [TUl E] 
which uses exactly the same setting. The difference arises at the following points: 

• The main difference is that our approach is not limited to the uniform 
baryon density distribution and takes the finite-size effects into account. In 
the limit of the uniform distribution, our approach will reduce to that of 

Refs. [iniin]. 

• The electromagnetic effects are not considered in ^OJ although the 
existing framework allows for accommodating them without further modi- 
fication. Since in the description of the pionic atoms where our approach 
will be applied, the Coulomb interaction plays an important role, we include 
the electromagnetic effects from the beginning. 

• Neglecting 4- and more-fermion interactions in Refs. fOl leads to sub- 
stantial technical simplifications. In our approach, we do not make this 
approximation in the perturbative expansion of the generating functional. 
Note that in ChPT there is no contribution from the 4-fermion interactions 
to the pion nucleus optical potential at 0{p^) (see below). For the treat- 
ment of the multinucleon interactions in the kinematical region where NN 
pairs are close to threshold, in addition, some kind of non-perturbative 
resummation might be necessary due to the large NN scattering length 

The formulation which was considered in Refs. jHl El E] does not use external 
sources. For this reason, in the uniform-density limit it is equivalent to our ap- 
proach only on the mass shell, i.e. in the calculation of the physical observables. 
Moreover, since Refs. [01 El Ej are based on the standard HBChPT description 
of pion-nucleon interactions the above equivalence holds up to the contri- 
butions from the non-standard counting regime (see below) ^. Finally note, that 
our approach which is based on constructing the S'-operator in the perturbation 
theory and then sandwiching it by the state vectors describing particles in the 
medium, is formally similar to the one used in Ref. fH] ^ again, in the case of 
the uniform medium. 

^Throughout the paper we use the terminology of Ref. for different momentum regimes 
in the Feynman integrals 
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We start with setting up the shorthand notation for the external sources in 
the generating functional 

] = {s,p,v,a}, Jo = {-M„ 0,0,0}, (2.1) 

where s,p,v,a stand for scalar, pseudoscalar, vector and axial- vector external 
sources, and Aig = diag (m„,md) denotes the quark mass matrix. The Green 
functions are obtained from the generating functional by differentiating with re- 
spect to the pertinent external sources, and setting j = jo at the end. 

The Lagrangian which we shall use to describe the system, is the standard 
Lagrangian of ChPT with pions, nucleons, photons and external sources in the 
vacuum: 

C = C, + Cf^ + £f ) + £f ) + Cf^'^ + ... + CP + C^^') + + . . . , (2.2) 

where the dots stand for the higher-order terms in the chiral expansion, for the 
terms of order and higher that are consistently neglected from now on, and 
for the terms in the Lagrangian containing four and more nucleon fields. For the 
convenience of the reader we collect the known expressions for the lowest-order 
chiral Lagrangians in appendix 1X1 

The asymptotic {in- and out-) states of the theory include both "elementary" 
particles (nucleons/antinucleons, pions, photons) as well as bound systems of 
nucleons (stable nuclei). On the other hand, bound systems of pions with nuclei 
- pionic atoms - are not present in the set of the asymptotic states, because in 
the present study their width can not be neglected. Instead, the pionic atoms 
will show up as (complex) poles on the second Riemann sheet in the pertinent 
Green functions. In this section, however, we consider the states containing 
only "elementary" particles. For example, the state |n; n; m; k; m(ont)) with n 
nucleons, n antinucleons, m pions and k photons, is given by 

IPiSi, ■ ■ ■ p„s„; piSi, ■ ■ ■ pnSn] Qi, " " " Qm; li^^i, " " " Ifc^^fc; in{out)) . (2.3) 

Here, Sj (sj) stands for the spin projection of the nucleon (antinucleon), and Si 
for the photon polarization vector. To ease notations, we suppress the indices 
that distinguish the different species of the nucleons (antinucleons) and pions. 

The generic off-shell transition process X — > F in the presence of any number 
of elementary fermions and bosons in in- and owt-states, with the sets X, Y 
both containing only bosons (pions and photons), is described in ChPT by the 
following matrix element 

(n'; n'; m'; k'; out\Toi{yi) ■ ■ ■ o^{y^)oi{xi) ■ ■ ■Op(xp)|n; fi; m; k; m) . (2.4) 

Here two sets of bilinear quark currents with pertinent quantum numbers (pseu- 
doscalar current for pions, vector current for photons): Oi(xi) ■ ■ ■ Op(xp) and 
Oi(yi) ■ ■ ■ Ouiiyuj) describe the bosons present in X and y, respectively. 



6 



The expression fj2.4|) is given in the Heisenberg picture. We further wish to 
rewrite this expression in the interaction picture. This can be done in a standard 
manner, using the evolution operator (see, e.g. [25). The resuh is given by 



(n';n';m'; k'\Toi{yi) ■ ■ ■ o^{y^)oi{xi) ■ ■ ■Op(xp) S\n; n;m; k) , (2.5) 

where S denotes the conventional S-operator, oi{yi) ■ ■ ■ o^{yi^)oi{xi) ■ ■ ■ Op{xp) 
stand for the currents in the interaction picture, and the free-particle states 
\n]n]m]k) are the eigenstates of the unperturbed Hamiltonian. Using Wick's 
theorem, the T-product in Eq. ()2.5|) can be expressed in terms of the normal 
products of free nucleon, pion and photon fields 

*,(x) = , (e-^^"fe,(k, s)u,ik, s) + e-^-"4(k, s^k, s)) 

sJ (27r)32yM2 + k 

$a(x) = I ^, ^ (a.(k)e-^- + at(k)e^^-), 
J (27r)32wVnf+k2 




A,{x) = E/ (^;^(^A(k)5,(A)e-^'=- + 4(k).;(A)e^'=- (2.6) 

Here, i = p,n and a = ±,0 stand for different nucleon and pion species. In the 
subsequent expressions we shall suppress these indices, if this does not lead to the 
confusion. Furthermore, Mj(k, s), fj(k, s) denote conventional free Dirac spinors 
and efj_{\) is the photon polarization vector. The state vectors from Eq. ()2.5|) are 
obtained by repeated action of creation operators on the perturbative vacuum 

n n m k 

In the canonical formalism, the free fields introduced above coincide with 
the interpolating field operators at time t = 0. At an arbitrary t, the relation 
between free and interpolating fields is given by the evolution operator. However, 
in order to avoid imposing boundary conditions on the interpolating pion field 
which is not uniquely defined, we reformulate the framework in terms of the 
functional integral depending on the external sources j. To this end, we consider 
the generating functional in the theory given by Lagrangian fl2.2j) 

Z{j\ri, V) = J VUVi)Vi:VA^ expji J {C + f]tp + ipri)d^x^ , (2.8) 
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where L = C{U, ip, ip, A^; j). Here U stands for the pion interpolating field matrix, 
A^ is the photon field, ijj = (^^^^ denotes the two-component interpolating 

nucleon field and * = #,„ its coujugate, , = ("') .s the nueleon external source 

and f] its conjugate. Note that the meson Green functions which are obtained by 
differentiating the quantity Z with respect to the bosonic external sources j, do 
not depend on the parameterization of the pion interpolating field U, the latter 
being merely an integration variable in the functional integral. Thus, the "off- 
shell ambiguity" mentioned in the introduction, never arises in this formulation. 

At the first step, we expand the generating functional ()2.8|) . in powers of 
nueleon sources r], fj 



Z{j\r^,fi) = Z^'\j) + J dUSjr]{x)Z^'\j\x,yUy) 
1 f 

+ ^ J d^xid%d\/idSj2f]{xi)f]{x2)Z^'^\j\xi,X2,yi,y2)viyi)viy2) ^ -(2.9) 

(Note that, since the interactions with the external sources j preserve the baryon 
number, we always have equal number of ?7's and f/'s in this expansion.) 
We further define the operator S{j) through the following construction: 



S{j) = 5(°)(j) + J d'xdSj : ^>{x)S^'\j\x,y)^>{y) : 

+ ]^ j d^x^d%dSjA ■■ ^(Xi)^^(x2)5(2)(j|xi,x2,|/i,|/2)*(l/i)^(y2) : + ■■■ ,(2.10) 
where 

S'^'\j\x,y) = zj,'" M) Z^'\3\x,y) {-I ^- M) z-^"\ (2.11) 

and so forth. Here M = diag (Mp, M„) and zn = diag {zp, z^) are, respectively, 
the 2x2 nueleon mass matrix and the nueleon wave function renormalization 
matrix. Furthermore, ": (■ ■ ■) :" denotes the normal ordering of the operators. 
The operator S{j) has the meaning of the S'-matrix defined on the subspace of 
the state vectors containing only fermions, in the presence of the external bosonic 
sources j. 

In order to obtain the matrix element (j'2.5p . we further differentiate the op- 
erator S{j) with respect to the pertinent bosonic sources, and expand the result 
in powers of j. This result is symbolically written as 

- [yi- ■ -Vu^xi- ■ -Xp) 



^jiiyi) ■ ■ ■ ^juj{yu>)^ji{xi) ■ ■■6jp{xp) 
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+ J d^zj^{z)Si^'^^'''{yi ■ ■ -yuj-^xi- ■ ■ Xp] z) 

+ \ j d'^zid^Z2ja{zi)jsiz2)S2^'^^'''\yi ■■■yuj]Xi--- Xp\ zi, ^2) H ,(2.12) 

where the index of the external source j corresponds to the quantum numbers of 
the quark current: jo- = v, a. Further, in a complete analogy with Eqs. ()2.9|1 . 
()2.1()|1 and (j2.11|l . we construct the operator 

S^^'^\yi ■■■yuj\xi---Xp) = Si^'^\yi ■ ■ ■yu.-.xi- ■ ■ Xp) 

+ j d^z: <l>„(z) : S^'^^'°'{yi---y^\Xi---Xp;z) 

+ j d^z: A^{z) : ^r'''^'''(l/i ■ ■ - y.; xi ■ ■ - x,; ^) + ■ ■ ■ , (2.13) 

where, in analogy with Eq. ()2.11|) . the coefficient S"^^'^^ " in ()2.13|) are obtained 
from S^^'-^^'" appearing in Eq. p.l2|l : e.g., S'^^'^^°'{yi ■ ■ ■ y^j] Xi ■ ■ ■ Xp] z) is ob- 
tained from iS^^'^^^d/i ■ ■ ■ i/i^; Xi ■ ■ ■ Xp; 2) by removing the pion pole in the Fourier 
transform with respect to the variable z, and multiplying the result by the wave- 
function renormalization factor ^Jz^. Integrating with the solution of the free- 
field equation $0(2) finally puts the 4-momentum, conjugate to the variable 2, on 
the mass shell. Other coefficients in the expansion ()2.13p are obtained similarly. 
We do not repeat their (evident) explicit expressions here. At the end, the oper- 
ator S^'-^^{y\ ■ ■ ■ yoj'-i Xi ■ ■ ■ Xp) is expressed in terms of the normal products of all 
free fields : 'I' ■ ■ ■ ■ ■ ■ ■ ■ ■ \I' : (we recall that S^'^''^^'" (yi ■ ■ ■ y^j] Xi ■ ■ ■ Xp] ■ ■ ■) 
are given by infinite series over normal products of fermion fields, cf. with Eq. 

Finally, the matrix element of the T-product of the currents in the Heisenberg 
picture (|2.4j) can be rewritten in the form of a matrix element of the operator 
S{Y,^} between the free-particle states 

{n';n'; m'; k'\S^^'^\yi ■■■yu>]Xi-- ■Xp)\n] n;m, k) . (2.14) 

This is the result we were looking for. The following remarks are in order: 

i) The construction presented in this section is the conventional LSZ formalism 
for the calculation of the matrix elements ()2.4p in the language of functional 
integrals (see e.g. [22]). Note that here one does not impose boundary 
conditions on the interpolating fields at t = 0. This is crucial for ensuring 
chiral symmetry at every step, and for circumventing the problem of the 
"off-shell ambiguity" . 

ii) The coefficient functions of the operators S{j) and S*^^'^^ are ultraviolet- 
finite, since the effective Lagrangian (j2.2p includes the set of counterterms 
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that make all S'-matrix elements finite. On the other hand, Z{j\ri,f]), in 
general, can not be made ultraviolet-finite with chirally-invariant countert- 
erms alone [25] . 

iii) In the presence of photons, the coefficient functions in S{j) and S^^'^^ are 
infrared- divergent. We shall regularize the infrared (as well as the ultravi- 
olet) divergences by using the dimensional regularization. The mass-shell 
limit is then understood in space-time dimensions different from four. 

iv) The coefficient functions in S{j) and S^"^'^^ are, in general, gauge-depend- 
ent, since they involve off-shell legs. Only those quantities that correspond 
to physically observable processes should be gauge- independent. 

v) In the calculation of the coefficient functions, some prescription is implied 
that preserves chiral power counting in the presence of nucleons. For ex- 
ample, one could have chosen heavy baryon ChPT ^3], or the infrared 
regularization method j2ni- As the most general prescription, we propose 
to use the so-called threshold expansion of the relativistic Feynman inte- 
grals [221 which, by construction, allows to keep the contributions from all 
soft integration regimes. In particular, the threshold expansion will au- 
tomatically produce the contributions to the Feynman integrals not only 
from the "HBChPT regime", but from the "non-standard regime(s)" fl] 
as well. Below we assume that the threshold expansion is always applied 
to all Feynman diagrams, even if this is not explicitly stated. 

2.2 Scattering processes in the presence of the nucleus 

The asymptotic spectrum of the theory described by the Lagrangian ()2.2|) in- 
cludes one-particle state vectors representing the nucleus as a bound system of 
A nucleons: 

|ri; in) = \n; out) = \n) , (ri'|0) = (27r)3 2P^ 6%P'^ - Pq) ■ (2.15) 

At this stage |f2) can still be any highly excited nuclear state as well. The 
scattering process Q + X ^ Q' + Y is described by the following matrix element: 

(f2'|Toi(yi) ■ ■ ■ o^(2/^)oi(xi) ■ ■ ■ Op{x,)\n) 
= {n'\ToM ■ ■ ■o^(y^)oi(xi) ■ ■ ■Op(x,)S|fi) . (2.16) 

Here 

= u(-oo,o)|ri) , (^]'| = (i^'|u(o,+oo) , (2.17) 

where \J{t',t) stands for the conventional evolution operator (see, e.g. |21j). 



10 



The basis |n; n; m; k) of the eigenstates of the free Hamiltonian is complete. 
A general nuclear state \VL) can be written as a linear combination of the basis 
states 

1^) = H Cn;n;m;k\n; n; m; k) . (2.18) 

n;n;m;k 

From the comparison to the expression of the bound-state matrix elements in 
the Mandelstam formalism, one can obtain the relation to the Bethe-Salpeter 
amplitude for the bound state 

Cn;n;m;k = J d'^^ ■ ■ ■ (^^^aG n;n,m;k{^l " " " X^) (0| *(0, Xi) ■ ■ ■ *(0,XA)|^1) , (2.19) 

where ^'(t, x) denotes the nucleon field operator in the Heisenberg representation. 
To make the interpretation easier, we found it useful to give this relation in a 
form containing the equal-time amplitude, which in the non-relativistic limit 
reduces to the Schrodinger wave function of the nucleus. In the language that 
uses the external sources, the Bethe-Salpeter amplitude is defined via the residues 
of the pertinent Green functions, without referring to the matrix elements of 
the interpolating fields. Further, Gn;n;m;fc(xi ■ ■ - ^a) stands for the perturbative 
kernel, which can, in principle, be calculated in ChPT. However, in our approach 
one does not need the explicit knowledge of neither this kernel, nor the Bethe- 
Salpeter amplitude. 

The quantity ()2.16|) which we want to calculate is given in the form of a matrix 
element of the operator S^^^^ introduced in the previous section, between the 
nuclear state vectors ()2.18j) 

(ll'lroi(yi) ■ ■ ■ o^(t/^)oi(xi) ■ ■ ■ Op{x,)\Q.) = (f]'|5^^^>(yi ■ ■ ■ y., xi ■ ■ ■ x,)\n) . 

(2.20) 

The right-hand side of this equation has the form J2iSi{^'\ : Pi : where 
Si denote the coefficient functions, and : Pj : stand for the normal products of 
the free fields : ^ ■ ■ ■ ■ ■ ■ -^/x ■ ■ ■ ^ '■■ The main idea of ChPT® consists in the 
following: one calculates the coefficient functions Si in standard ChPT, whereas 
the matrix elements of the normal products of the free fields are parametrized 
in terms of a few characteristic functions describing nuclear dynamics. This is 
the place where the empirical input enters in our calculational scheme^. Our 
strategy consists in relating these matrix elements to a limited number of nuclear 
observables, and predicting other observables where the same matrix elements 
occur. 

^As it stands, the approach is based on a well-defined separation of the perturbative and 
non-perturbative contributions in the matrix elements of currents. Using the same arguments 
one could, in principle, reshuffle particular contributions between these two groups, and design 
a different scheme. Here, we shall not explore this possibility any further. 
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If the approach is used to describe a large nucleus {A — > oo), in its ground 
state, the following assumption can be justified: 

Assumption 1. Matrix elements of operators containing at least one free pion 
or photon field, vanish between the nuclear states. We interpret this property in 
the following way: \VL) describes the ground state of the nucleus as a system of 
nucleons only aX A ^ oo. 

The above assumption enables one to rewrite the matrix element ()2.16|) in 
terms of the operator S{i) defined by Eq. ()2.10|) . as 



Sjiiyi) ■ ■ ■Sjc,iyc,)6ji{xi) ■ ■ ■ 6jp{xp) 



(2.21) 



J=JO 



Assumption 2. In order to simplify the present calculations, we assume that 
the nucleus is in the spin-saturated state. This means that the nucleus is not 
characterized by internal quantum numbers. The latter assumption, which can 
obviously be relaxed on a later stage, in combination with the static limit {A — > 
oo), helps to reduce significantly - on the basis of symmetry considerations alone 
- the amount of empirical input needed on the nuclear side. 

Assumption 3. Finally, we assume that the nuclear background does not 
affect the ultraviolet properties of the theory. In the case of the infinite nuclear 
medium, this is an exact statement, since all diagrams with medium insertions 
contain a cutoff at Fermi-momentum kp. The generalization to finite systems 
amounts to require that Fourier transforms of distributions involving \Q) drop 
sufficiently rapidly at large momentum. Put differently, we assume that in the 
expansion J2iSi{^l'\ '■ Pi '■ \^) each matrix element is finite, and there are no 
cancellations of divergences between various terms (we recall that 5*^ are finite 
after renormalization). 

It should be emphasized that, with the above assumptions, ChPT® can be 
pursued to any given chiral order. The consistency of the approach, and the 
validity of the above assumptions, should then be verified a posteriori order by 
order in the chiral expansion. 



2.3 Pion-nucleus bound state 

The primary aim of the present paper is to study states of negatively charged 
pions bound to heavy nuclei. To this end, it is useful to consider the two-point 
function of the pseudoscalar densities, sandwiched by the nuclear states 



, (2.22) 

j=jo 



where F is the pion decay constant in the chiral limit, B is related to the quark 
condensate {F^B = —limm^ ma^oiQl))- We follow the Condon-Shortley phase 
convention and express the isovector (pseudoscalar) source field as p^{x) = 
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(2)^-'^/^(=l=Pi(x) + ip2{x)), p^{x) = P3{x). The normalization in Eq. ()2.22|) is 
defined so that in lowest order of the chiral expansion, the residue in the pole of 
this two-point function is normalized to unity. 

Using translational invariance, we remove an overall (5-function, corresponding 
to the conservation of the total 3-momentum, from the quantity Dq/qIq', q). In 
the remaining matrix element one may pass to the static limit Pq = Pq = 
Mq{1, 0). In the calculations, it is convenient to use the Hilbert space where the 
static nucleus is described by the state vector \Q) normalized to unity ^. Using 
the new notations, the right-hand side of Eq. ()2.22j) at A — > oo can be written as 

Dt,,a(q', q) = (i'^fS^iK + q' - Pli - ^W, q) , 
DW,q) = /AdV"'— 

For a static nuclear background, energy conservation implies 

Diq', q) = 27r5(g°' - q', q) , E = gO' = . (2.24) 

The pion self-energy is defined from the inverse of the operator D as follows: 

U{E; q', q) = i27ifS%q: - ci)D^\E; q^) - D~\E; q', q) , (2.25) 

where Dq^(E; q^) = E'^ — — is the inverse of the free pion propagator, and 
= mT^+ stands for the physical mass of the charged pion in the vacuum. 
Our strategy for finding the energy levels and widths of the pionic bound state 

implies the following steps: 

1. Calculate the non-local pion self-energy operator Il{E; q', q) in the system- 
atic ChPT expansion. 

2. Find the pole position in the two-point function Z)(£';q', q) with a given 
self-energy operator Il{E; q', q), using non-perturbative methods. The real 
and imaginary parts of the pole position are related to the real energy and 
decay width of this state. 

3. At the last step, for a given self-energy operator calculated perturbatively 
at a fixed order in ChPT, study the chiral expansion of the energy and the 
decay width. This enables one to control the convergence in the bound- 
state characteristics at all stages of calculations, possibly eliminating (po- 
tentially large) higher-order contributions that arise due to the use of the 
non-perturbative techniques. 

•^This state vector, which describes the nucleus in the static Umit. is in fact defined by Eq. 
(|2.23|) . The more detailed treatment of the static limit is considered below, in section 



, (2.23) 
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It should be understood that identifying the (complex) bound-state energy 
with the pole position of the pion two-point function in the nuclear background 
(12 .2411 implies additional assumptions: 

Assumption 4- In Eq- fl2.22|l the nucleus is on shell in its ground state, 
whereas in the pion-nucleus bound state both pion and nucleus are off shell. The 
rigorous definition of the bound-state energy involves finding the pole position of 
the Green function of 2A fermion and two pion fields. We expect, however, that 
for a very large A one may safely put the nucleus on mass shell. 

Assumption 5. If the nucleus is put on shell, this means that the effect coming 
from nuclear excited states is neglected. In order to demonstrate this, note that 
the initial-/final- state interactions in the A ingoing/outgoing nucleon legs of any 
Green function create a tower of radial-excited nuclear bound states. Putting 
the nuclear momentum on the mass shell eliminates all but the ground state 
contribution to the Green function and thus, from the mathematical viewpoint, 
amounts to "freezing" the nucleus in its ground state during the whole interaction 
with the projectile pHOlj . 

By construction, the energy of the bound state obtained with this procedure, 
does not have any spurious dependence on the choice of parameterization for the 
pion field U. If, instead of the quantity D given by Eq. ()2.22|1 . one considers 
the two-point function of the interpolating pion field, such a dependence arises 
in this two-point function itself ^21 EI- As one would expect, this dependence 
indeed disappears in the pole position up to a given order in the chiral expan- 
sion When solving the bound-state equation numerically for the non-local 
self-energy jT2], it is in general difficult to demonstrate that the dependence on 
the parameterization of the pion field is of higher chiral order. It is therefore 
advantageous to work in a framework which avoids this dependence ah initio [of 
course, the pole positions in the two-point functions of the pseudoscalar densities, 
and of the pion fields on the other hand, should coincide order by order in the 
chiral expansion.]. 

We finally note that the two-point function is gauge- dependent. In order to 
demonstrate that the approach is consistent at a given chiral order, one has to 
explicitly show that the pole position is invariant up to this order. In this paper, 
we check the gauge invariance up to and including 0{p^) (see section 121). At 
present, we are however not aware of the general proof of the gauge invariance 
that extends to all orders within our approach. If it turns out at higher orders, 
that the gauge invariance of the bound-state energy is lost, this would signal about 
the inadequacy of the approximations which were used in the formulation of the 
present framework. From the phenomenological point of view, the issue is totally 
irrelevant, since it concerns only the higher-order electromagnetic corrections to 
the pion-nuclear optical potential. 
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2.4 Free Fermi- gas 

In order to examine the limit of a homogeneous nuclear medium, it is instructive 
to start with the simple case of a Fermi gas. The ground state of this system 
has all levels filled with nucleons up to the Fermi- momentum kp- To simplify 
notations, we assume here in addition that isospin symmetry is not broken, switch 
off the electromagnetic interactions, and suppress the index labeling different 
nucleon species p and n. The ground state in this approximation is 

\n)^\^l,)=U n ^^(l.,5.)|0), (2.26) 

\i,Si=±l/2 

where |0) denotes the perturbative vacuum, and A/" is a normalization constant. 
Furthermore, we again introduce the "round-bracket" notation, with (f2o|f2o) = 1 
[see section 1221 for further details.]. 

Using the canonical commutation relations, one can easily show that 

(fio|&^(li,si)Kl2,S2)|l^o) = {2nf2ll5\h - h)5s,s,e{kF - |Ii|) , 

(no|&^(li, s,)h\h, S2)6(l3, s^)h{h, s^)\n,) 

= (27r)3 2/0 5\h - h)6s,s, e{kp -\h\) {27rf 2/° 5^{h - h)5s,s3 dikp - \h\) 

- [perm. 1^2], (2.27) 

and so forth. The free fermion (nucleon) propagator in momentum space, with 
Qq in the background, is given by 

= I J dW^^ (0|T^(x)^^(0)|0) + I J d'^e'P '' {Qo\ : ^(x)^(O) : \Qo) 

= (M^+ ^) { j^2_\2_,Q - ^""'^(P' - Ml)e{py{kF - IpD) , (2.28) 

where M,j, stands for the mass of the fermion field. One immediately recognizes 
that the quantity S^i^ [p) is the relativistic in-medium free fermion propagator. 

In the calculation of the Green functions of the external sources j, one needs 
to evaluate only the closed fermion lines at non-zero fermion density (see Fig. 
Let us first consider the closed loop with two vertices ri^2, where these stand for 
arbitrary matrices in the space of Dirac and flavor indices, describing the coupling 
of the external sources j to the fermion bilinears. Applying Wick's theorem and 
using Eq. ()2.27p for the normal products, we get 

= -Tr [Snoix - y)T2Sno{y - x)T,] + Tr [^^0(0)12] Tr [SnM^i] ■ (2-29) 
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Figure 1: Closed fermion loop with any number of vertices Fi, F2, ■ ■ ■. Double 
lines denote the external sources j attached to the loop. 

The above formula states that the closed fermion loop with two vertices in the 
background of is calculated replacing free nucleon propagators by the in- 
medium fermion propagators (x). The same rule holds for a loop with any 
number of vertices, or any number of factorized loops. Indeed, we may redefine 
the normal product with respect to the new vacuum \Qo). Then, according to 
Eq. ()2.28|1 . the fermion propagator is just S'qo(x). Since the combinatorial factors, 
produced by the contractions in the Wick's theorem, do not depend on the explicit 
form of the propagator, one finds that all differences between closed fermion loops 
calculated in the vacuum and in the Qq background reduce to the replacement 
of the nucleon propagator in the vacuum by Sqq{x). One concludes that, for a 
uniform fermion density distribution Qq, the approach described in the present 
paper reduces to the standard in-medium ChPT. 

2.5 Chiral power counting 

Another crucial ingredient for the approach proposed here is the systematic book- 
keeping based on chiral symmetry. On first sight, the formalism allows too much 
freedom in the nuclear matrix elements, which have to be evaluated for arbi- 
trary space-time arguments of the fields We shall now demonstrate that 
the chiral power counting, in general, enables one to substantially reduce this 
freedom. 

Even for a non-uniform distribution of baryon number, it is still useful to 
think in terms of the "local Fermi momentum kp" which in this case is a function 
of the space coordinate. It is common and convenient to count kp at 0{p) [the 
actual numerical value for this quantity in heavy nuclei, kp ^ Sm^r, supports this 
conjecture]. Consequently, the kinetic energies of the nucleons can be counted 
as 0{p'^). The fermion field itself is 0(p^/^), ~ 0(1), and the counting in 
the coefficient functions S**-"^ is standard. Furthermore, in order to be conform 
with the chiral power counting for the pion self-energy operator (j2.25|) in the 
uniform medium, we declare this quantity to be calculated at 0{p''), if its actual 
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power according to the above counting is k — 3 (we remind the reader that, 
for the uniform distribution, the self-energy operator reduces to Il{E; q', q) 
(27r)^(5^(q' — q) Il{E] q^), and k thus corresponds to the chiral power of Il{E] q^), 
whereas three additional powers are "eaten" by the (5-function.). 

The chiral power counting described here implies a systematic expansion of 
the matrix elements in Taylor series, so that finally one ends up with the matrix 
elements with fields ^'(a:) and derivatives thereof, taken at coinciding space-time 
arguments. Such matrix elements can in turn be written in terms of local distri- 
butions. This makes their construction transparent and reliable. 

Below we shall schematically explain this method for the case of the normal 
product of two fermion fields. The corresponding contribution to the vr-nucleus 
scattering matrix element is given by 

= J d^d\jD^^^{q',q\x,y) : t,(x)vl>,(y) : , (2.30) 

where 

U,, (g,g|x,yj - ^^^^^^ j d^dze Sp+{u)5p~{z) ,=,o ' ^ ^ ^ 

and a = («,z), b = (/3,j) stand for multiindices with a,P and i,j denoting 
Lorentz and isospin indices, respectively. Furthermore, the threshold expansion 
is implied only for the quantity S**-^-* which is calculated in the vacuum. There is no 
need to perform the threshold expansion in the matrix element {Q\ : ^b{x)'^a{y) '■ 
since, according to our picture of the nucleus, all 3- momenta involved receive 
a cutoff at kp. 

Let us first consider the case of the standard HBChPT counting. In this 
case, both energy and 3-momentum of the pions in the Feynman diagram count 
like 0{p). In order to express the nuclear matrix elements in terms of local 
distributions, we introduce the center-of-mass and the relative coordinates x = 
R + s/2 and y = R — s/2. Since the nucleus is static, the matrix element does 
not depend on the time component of R. Now, we consider the dependence 
of this quantity on s°. Factoring out the nucleon (proton or neutron) mass Mq,, 
which is counted at 0(1), we can expand the rest in the Taylor series of 

{n\ : ^,{x)^a{y) ■■ \n) = e*^^*^° {n\ n,{o, x)iv„(o, y) : \n) 

_ e^M^.o ^ . ^^^Q^ - ^^^Q^ y) . |(^) + . . . ^ (2.32) 

where Na{x) = e*^^*^° "^aix). Every time derivative acting on Na{x) produces 
the kinetic energy of the nucleon, and has therefore to be counted as a quantity 
of order p"^. On the other hand, every power of s° corresponds to the differentia- 
tion with respect to the energy variable in the Fourier transform of S^^\ Since, 
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according to the standard HBChPT counting, differentiation of every fermion 
propagator reduces the chiral power by 1, one may count s° at 0{p"^). We fi- 
nally see that every term in the expansion ()2.32|1 is 0{p) as compared to the 
previous term, so, at a given order in the chiral expansion, one may truncate the 
series, and arrives at a result that is local in the relative time variable s°. 

It is clear that this argumentation does not work for the spatial nonlocality: 
here, every space derivative in the matrix element counts as 0{p), and spatial 
components ~ 0{p~^), so one is not allowed to truncate the Taylor series. 
One may, however, expand the fermion propagators contained in S^^\ In the 
HBChPT regime, this expansion reads 

- A±£ + J_ + il±^] + o(^) (2 33) 

where p = Mq,v + /, the unit vector v = (1, 0), and as usual ^ = p^'~i'^. From this 
expansion one observes that i) each term with an additional power of M^^ is 0{p) 
suppressed with respect to the previous term, so the series can be truncated at a 
given chiral order, and ii) each term is a polynomial in 3-momenta. According to 
the above properties, at a given order in the chiral expansion, S''-^^ in coordinate 
space is a finite sum of spatial (5-functions and derivatives thereof. Substituting 
into Eq. 1)2.301) and integrating by parts, one finally obtains a finite sum containing 
matrix elements of the fermion bilinears and derivatives thereof, taken at the same 
space-time point R. In addition, we recall that the field ^ obeys the free Dirac 
equation, so all time derivatives can be eliminated at the end. 

The expression which one obtains after performing all Taylor expansions, is 
much simpler than the original one. Whereas the matrix element in Eq. ()2.30p in- 
volves a "dynamical" density matrix (with non-coinciding time variables), which 
is strongly dependent on the description of the structure of the nucleus, after the 
Taylor expansion one gets terms of the type 

(fi| : f b(0, K)K,--- i *a(0, R) : |1^) , n = 0, 1, • • • , (2.34) 

which have a simpler interpretation. For example, using the Fierz transformation, 
the term without derivatives can be related to the matrix elements of local scalar, 
vector, ... fermion bilinears {VL\ : ^Ti"^ : \Q) with Fj = 1,7^, To lowest 
order in the chiral expansion, these matrix elements coincide with the scalar, 
vector . . . formfactors of the nucleus, which can be parametrized, e.g., in terms of 
the pertinent radii, without any more detailed explicit knowledge of the nuclear 
structure. 

Next, we briefly consider the case of the non-standard counting. As it will 
be demonstrated below by explicit calculations, the contributions to the self- 
energy, corresponding to the non-standard counting regime, do not appear up to 
and including 0{p^) in the pion two-point function in the vicinity of the mass 
shell. For the general pion kinematics this is not always the case. In particular. 
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it has been shown [2H1 that, for a large class of the collective phenomena, the 
elementary excitations can be described by the non-relativistic dispersion law 
E = 7p^. Moreover, it has been argued, that the same dispersion law applies 
in the Goldstone boson condensed phase of the nuclear ground state (see e.g. 

where this issue is considered in detail). So, in order to be able to use 
ChPT® to study the phenomenon of the pion (and kaon) condensation in general 
[21 EOl ISl E21 ! it is necessary to modify counting rules for the external momenta 
and to assign ~ to the energy, while 3-momentum counting remains standard 
[TT] . At the next step, one has to look for the pole in the two-point function 
with the vanishing energy E, which signals that the condensation has indeed 
occurred. One sees (section 4.2 of Ref. jTJ) that in the asymmetric matter 
where the density exceeds the critical value given by Eq. (4.10) of this paper, 
the inverse pion propagator indeed develops a pole at i?, — > and D = p'^/E 
finite. This means that for the densities that obey the condition (4.10), the non- 
standard counting may appear already at the tree level. Note as well, that Eq. 
(4.10) implies that the density is counted at 0{p^) rather that at 0{p'^). Another 
example, where the non-standard counting appears at the leading order, is the 
in-medium tttt scattering amplitude In this paper, we do not consider the 
non-standard regime in detail since in the context of pionic atoms the contribution 
from this regime arises at the one-loop level. In this case, our approach should 
be adapted correspondingly. 

To summarize, chiral power counting helps to systematically organize the 
perturbation expansion in the background of a static nucleus, in terms of local 
distributions. This construction, however, heavily relies on the assumption that 
only the momentum regime relevant in HBChPT, gives rise to non-zero contri- 
butions in the coefficient functions S^'^\ This conjecture may or may not hold 
- in the latter case one has the non-standard counting which has to be dealt 
with separately. Non-local quantities, describing more detailed aspects of nuclear 
structure, may become necessary here. 

In the section |21 we present the explicit calculation of the pion self-energy 
operator at 0{p^). As already mentioned above, we demonstrate that only the 
HBChPT momentum regime is relevant in Feynman integrals at this order, so 
the answer can be given entirely in terms of local distributions. In this paper we 
do not attempt to analyze higher orders in the chiral expansion. 

2.6 The static limit 

In this section we examine in some detail the reasoning behind the picture of a 
static nucleus, and explain the meaning of the "round-bracket" notation which 
was introduced in Eq. (j2.23|l . As was stated, the procedure involves two steps 
when evaluating the matrix elements (j2.22j) in the background of a (heavy) nu- 
cleus: first, removing a momentum conserving delta function in the CM frame 
[which coincides with the nuclear rest frame for an infinitely heavy nucleus] , and 
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secondly, passing to the static nucleus limit in terms of the velocities of the initial 
and final nuclei. 

For the sake of illustration, let us first consider the term in the two-point 
function ()2.22|) . corresponding to the in- vacuum self-energy of the pion: 



3=30 

(2.35) 



with = 1. 

Next, consider the term with one medium insertion. Using translational in- 
variance and taking the limit Mq — > oo, we obtain [cf Eq. (j2.30|) ] in the CM 
frame: 

J d'xdSjDi'J{q',q\x,y){n{P')\ : f,(x)vl>,(y) : \n{P)) 

= (27r)353(p' + q'_p_q) J dU^oi^J {q' , q\x,y){n\ : ^,{x)<i/a{y) : 1^) , 

(2.36) 

where 

(fi| : ^t,ix)^aiy) H^) = / ^ e^"^-^ /,,(Q, X - y) , 

/.a(Q,x - y) = (fi(-L - |)| : t,(0, ^)*a(o, : M-L + |)) , 

(2.37) 

and L = | (q' + q) [Note that in the static limit oo the bracket matrix 

element in Eq. fl2.37|) should not depend on the momentum L, but only on the 
momentum transfer Q since, in this limit, the round-bracket matrix element in 
the same equation should depend on the vectors x, y alone.]. From these examples 
we arrive at the following interpretation. The round-bracket state \Q) denotes 
the static nucleus fixed at a given point of space (at the origin). As mentioned 
already, this state is assumed to be spin-saturated for simplicity. The expectation 
value of any operator in the state \Q) is translation-invariant in time (but not 
in space). Furthermore, this expectation value (round-bracket matrix element) 
is equal to the Fourier transform of a corresponding formfactor evaluated in the 
Fock space. The generalization for the higher-order terms is straightforward. 

Next, we consider in more detail the role of the restriction to spin-saturated 
states Q. Let us start with the matrix element of the operator : ^(x)7^\l/(x) :. 
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In order to pass easily to the static limit, it is advantageous to work in terms of 
velocities = P^/Mq, v'^ = P'^/Mq instead of momenta P^, P'^. The above 
current is conserved since is a free field. Then, from the assumption that the 
bound state of the nucleus does not depend on the internal quantum numbers, 
it follows that the Fock space matrix element of this current is described by a 
single scalar formfactor: 

(n(P')l : ^^(0)7,^(0) : =^-{v' + v\F{t) , t = {P' - Pf , (2.38) 

where F{t) is the form factor of the charge distribution related to the current 
In the static limit, when w^, v' ^ (1, 0) and t -(P'-P)^ = -(q'-q)^, 
according to Eq. ()2.37|) . this matrix element reduces to 

(^]| : ^(x)7,^(x) : = 9,0 J e^^^ F(-Q^) . (2.39) 

Other fermion bilinears can be considered analogously. For example, symmetry 
considerations imply that the matrix element of the operator : ^{x)a^u'^{x) : 
can involve only the Lorentz-structure v'^^v^ — v'^v^. Then, in the static limit, 
the round-bracket matrix element should vanish: 

(fi| : t(x)a^,^(x) : |fi) = 0. (2.40) 

To summarize, the use of equations of motion and symmetries, together with 
the simplifying choice of a spin-saturated background nucleus f2, greatly reduces 
the number of independent round-bracket matrix elements which serve as an 
empirical input in the construction of the self-energy operator. 

Finally, we mention that for the free Fermi-gas the equality (f2o|^o) = 1 and 
further, the equations (|2.27p in the rigorous sense are understood as follows 

(fio|&^(li,si)Kl2,S2)|fio) = {2Txf2ll5\h-h)e{kF-\\,\){n^\no) ,(2.41) 
and so forth. 

3 Pion self-energy at 0{p^) 

In this section we demonstrate the general rules formulated in previous sections, 
by presenting the detailed calculation of the pion self-energy operator at O(p^) 
in ChPT in the presence of the finite nucleus. This is an instructive exercise, 
even though the calculation up to this order yields just the leading terms of the 
pion-nuclear optical potential, those linear in proton and neutron densities. 
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3.1 Pion self-energy in the vacuum 

The two-point function of the pseudoscalar densities in the vacuum is given by 
[cf with Eqs. ^TWf and ((2221)] 



(3.1) 



Up to and including O(p^), only the diagrams shown in Fig. |21 contribute to the 
self-energy of the pseudoscalar densities in the vacuum. The calculations are 
most easily done in the so-called a-model parameterization 

(3.2) 



F2 

where 7^{x) denotes the interpolating pion field. We use this parameterization 
throughout this paper [of course, the two-point function of the pseudoscalar den- 
sities does not depend on the choice of a particular parameterization]. A straight- 
forward calculation with the use of the Lagrangian from appendix |X] gives 

D^'\^') = + ^^n(°)(g^)^^ + . . . , (3.3) 

where the dots stand for higher-order terms in the chiral expansion, and 

nW(g2) ^ (q^-mDA'-'^ + iq'-mlfB^'^q'), (3.4) 

where B^^\q^) is finite at q^ and d ^ A. The pertinent contribution to the 

self-energy is given by 

n(°)(E;q',q) = {2'Kf6\c{ - q)n(°)(E2 - q^) . (3.5) 

The expression for A^^^ takes the form 

= 7l(2'4-I^'"f)+^'3 + ^'|2(3-m,„ 

- ^(*:; + *:^-2*J-2*S) + ?*; + 8*a+0(rf-4), (3.6) 

where n denotes the scale of the dimensional regularization. Note that we tame 
both ultraviolet and infrared divergences in this regularization, and we attach 
the subscript "IR" to 

in order to distinguish ultraviolet and infrared divergences. Finally, we note that 
the explicit expression for 5(0) (p2) 

is never needed, because it does not contribute 
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Figure 2: Diagrams contributing to the two-point function of the pseudoscalar 
densities at O(p^) in the a-model parameterization. Double, dashed and wig- 
gle hues stand for the pseudoscalar densities, pions and photons, respectively. 
Crosses denote the vertices with O(p^) LECs fcj. 

to the residue of the two-point function at this order of the chiral expansion - for 
this reason we do not display it here. 

As one immediately sees from Eq. ()3.6p . the pion self-energy is both infrared- 
divergent and gauge- dependent off-shell, i.e. when ^ m^. Below we shall 
demonstrate that in the eigenvalue equation for determining the binding energy 
of the pion-nuclear bound state, where pion self-energy plays the role of the 
potential, such off-shell effects begin to contribute only at higher chiral order and 
can be neglected at the level we are working. 

3.2 One medium insertion 

According to Eq. (j2.10|) . the contribution with one medium insertion is expressed 
in terms of the coefficient function (j2.31|) where, up to and including O(p^), the 
quantity D^^^ is given by the diagrams in Fig. |21 These are: the one-photon 
exchange diagram, contributions from the anomalous magnetic moment of the 
nucleon, the Weinberg- Tomozawa vertex, contributions from the LECs Cj, /j, and 
the nucleon pole diagrams 

{q" - ml){q' - ml)Dil\q', q\x, y) = 5\x - y)S'^'~''y^[D\^:\q\ q) 

+ Dt:\q\ q) + D^r\q'. Q) + (q' , q)] + DZ\q\ y) , (3.8) 

where we use the physical value of the squared pion mass, m^, instead of its 0{p^) 
value, since the difference is of higher order in ChPT, which we neglect anyway. 
Further, 

Dt'\q\q) = e'{^,),J-il + T%{q' + q),D^'^{q'-q), 

Dt\q'.q) = -(^^.AW(^^^ + f r^) (q' + qUq-qTD^'^iq'-q), 
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\ / \ / \ / ^^^^ 
> > i > > > W> — > > 



Figure 3: Diagrams contributing to the quantity D^^'^ up to and including 0(p^): 
(a) one-photon exchange; (b) contributions from the anomalous magnetic moment 
of the nucleon; (c) Weinberg- Tomozawa vertex; (d) contributions from LECs 
Ci, fi] (e) nucleon pole diagrams. 

2 



{4l,S,{x - y) /75)/3«2 (1 - r^),.e^'^'"-^ • (3.9) 

In these formulae, D^^'^ and Sp^n stand for the free photon and baryon propagators, 
respectively: 

= / T^e-^'-" -. (3.10) 

^ ^ i (27r)4 Mi- / ^ ^ 

The corresponding contribution to the self-energy is 

27r<5(g°'-g°)n«(E;q',q) = j d'xdSjD\^Xq\q\x,y){n\ : ^,{x)^a{y) : |^^) (3.11) 

We proceed to investigate term by term in detail the contributions (|3.9p to the 
two-point function of the pseudoscalar densities. 

3.2.1 Coulomb potential and contribution from the anomalous mag- 
netic moment of the nucleon 

According to Eq. (j3.9j) . the one photon exchange contribution to the two-point 
function of the pseudoscalar densities is given by 

D^''\q\q) = ^^^27r%°'-g°)n(i^)(i?;q',q)^^, (3.12) 
q — mz. q^ — 
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where 

27r(5(g°' - g°)n(^^)(E; q', q) = J d^xe'^"' -"^'^ {q' + q),D^''{q - g) x 

X (l]|:f(x)7^1(l + r=^)vl>(x):|fi). (3.13) 

In order to demonstrate that the longitudinal piece depending on the gauge pa- 
rameter ^ does not contribute in Eq. (j3.13j) . we write this matrix element in the 
form 

{Q\ : f (x)7,^ (1 + r')^ix) : \Q) = ^,o^ (p°(x) + /(x)) , (3.14) 

with the isoscalar and isovector densities = Pp + Pn and p^ = Pp — Pn- The 
density in ()3.14|) is related to the Fourier-transform of the formfactor of the 
nuclear charge distribution, 

- (p°(x) +p3(x)) = p,(x) = J ^e''^^o+3[-l^] , (3.15) 

where 

(O'(P')I : ^(0)7^(1 + ^')^(0) : l^(P)) = I W + v) ,Fo^s[t] (3.16) 

with t={Pl^- PnY [cf Eqs. flTT^ ^^]. At the order in ChPT we are working, 
the formfactor Fo+3[t] is proportional to the electromagnetic formfactor of the 
nucleus Q 

The presence of S{p^ — q^) in Eq- ()3.12j) implies that the one-photon exchange 
contribution turns out to be effectively instantaneous. Eq. ()3.13p takes the form 

n(^^)(E; q', q) = -2E J d^^e-'^'^'-'^^ Wd^) , (3.17) 

where 

Vc{^)= fd'r^^ppir), (3.18) 
J |x — r| 

and a = e^/(47r) denotes the fine structure constant. Recall that, in the chiral 
counting, p^^'^^(r) ~ and e ~ p. It is straightforward to observe that the 
quantity n^^'^^ given by Eq. ()3.17|1 . counts as 0{p). According to our conventions 
from section 121 [the actual chiral power of Il{E; q', q) differs by 3 units from the 
chiral power assigned, in order to stay in conformity with the chiral counting in 
the vacuum], the Coulomb term contributes to the self-energy at O(p^). 

We finally note that the contribution from the anomalous magnetic moment 
of the nucleon - which otherwise would contribute at 0{p^) - vanishes since the 
matrix element of the operator : 'i'cr^t,\& : between the nuclear states disappears 
in the static limit [see Eq. ()2.40|) ]. Consequently, 

n("")(E;q',q) = 0. (3.19) 
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3.2.2 Contact contributions 



Contact contributions include: tlie Weinberg- Tomozawa term [0(p^)], and the 
LEC contributions [0(p^)] involving q and . From Eq. ()H.9|) one obtains 

jl(wn+(cf)^E; q', q) = I d'xe-'fo'-o'- |n„(£;; x) + q' • qn2(B; x)] . 
n„(E; X) = - J, A.) + ^ 4c....;.-2fe + C3)E'^ ^ .(^j ^ ^^^3(,) , 



n2(E;x) = ^a°(x), (3.20) 

with the isoscalar and isovector scalar densities 

a[°'3l(x) = : ^{x)\\,T^\^{x) ■ . (3.21) 

Note that the contribution proportional to C4, vanishes for the static nucleus. 
We would like to stress that, in order to be consistent with the chiral counting in 
the nuclear matrix elements and our description of the nucleus, one has to count 
the difference between a^^''^"^{y^ and p[''''^](x) at a higher chiral order than these 
quantities themselves. This difference, which involves 1 — 70, is determined by 
the overlap of the "small" components of the Dirac wave function of the nucleon. 
For a slowly moving nucleon inside the nucleus, this difference is suppressed by 
a factor as compared to the overlap of the "large" components. This can be 
directly seen for the case of the uniform density, where 

p[o,3] _ ^[0,3] _ 0{k%), 

whereas each term individually is of order k\. Below, we shall always use the 
counting 

^p[o,3i(,) ^ ^ - Oij?) , (3.22) 

and eliminate o"[°'^l(x) in favor of p[°'^l(x) in all expressions. 

Note also that one may relate the matrix elements 1)3.211) to the scalar form- 
factor of the nucleus 



J=Jo 



+ ••■, (3,23) 

i=io 



modulo higher-order terms in chiral expansion. At the order we are working, 
these higher-order terms can be neglected. 
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3.2.3 Nucleon loop in the presence of the nucleus 



The contribution coming from the diagrams of Fig. |Hfe, involves the matrix el- 
ement of two fermion fields taken at different space-time points. In order to 
evaluate this contribution, one has to use the expansion described in sections l2.5| 
12.61 along with the bookkeeping based on chiral symmetry. A straightforward, 
yet tedious calculation leads to: 

n(^)(E; q', q) = ^ / d'^e-'^'^'-'^-'iME; q', q|x) + U,{E- q', q|x)] + ■ ■ • (^.24) 
where 

n4(E;q',q|x) = ^p'(x), 

u fj^ ' I ^ (M„ - Mp)q^ ■ q 3 
n5(E;q,q|x) = p^(x) 

^ i^^-2q-.q , 3(q--q)^-qV ^ , 

Here Mjq = | (Mp -|- M„), and the dots in Eq. ()3.24p stand for the higher-order 
terms in chiral expansion. 

3.3 Two medium insertions 

The coefficient function D^'^\ which corresponds to the two medium insertions in 
the pion self-energy, is given by 



(25F)2 J 6p+{u)6p-{z) 



3=30 



Up to and including O(p^), there is a single contribution which stems from the 
pseudovector pion-nucleon Lagrangian. This contribution, which is diagrammat- 
ically shown in FiglH equals 



X 



+ crossed term . (3.27) 
The pertinent contribution to the self-energy is 



2vr5(g°' - q')U^'\E; q', q) = / d'xd\e^^' 



x-iq-y ^ 
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> ^ > > ^> 



Figure 4: Two medium insertions: the diagrams contributing to the quantity 
Eq. dSZZD- 

X {n\ : ^{x) /75(r^ + ir"^)"^ (x)^ (y) ^75(r^ - zr^)^(y) : \n) + crossed term. 

(3.28) 

It is convenient to work in the momentum representation for the free fields 
For the self-energy, one gets the following expression 



SIS2SSS4 



X Si) 4'l^u{\i + g', S3) u(l2, S2) #75^^(12 - g, S4) x 

X (fi|6t(liSi)6t(l2S2)6(l3S3)Kl4S4)|f^) + ' " " , (3.29) 

/ u. (Is) \ 

where m(Is) = ''^/i ( denotes the 8-component Dirac spinor of the nucleon, 

\Un{\s) J 

and the dots indicate terms with creation/annihilation operators for antinucle- 
ons, and the crossed term. Below, we shall demonstrate that, in analogy with the 
infinite nuclear medium the above contribution vanishes (terms with antin- 
ucleons can be considered similarly). The reason for this lies in the fact that 
we still imply a soft cutoff on the 3-momenta of the nucleons within the state VL 
[albeit with the position-dependent "local Fermi momentum" kp\. On the other 
hand, the argument of the (5-functions in Eq. ()3.29p does not vanish in the soft 
momentum region. For example, the solution for |Ii| is 



-A|q'|cos^± ^A2|q'|2cos2e- (M2g^," - A2)(g^' - Iq'^cos^^) 

|li| = 72 TTi^ —n ' (3-30) 

% — |q I cos"^ U 

where 9 is the angle between q' and li, and 

g'2_ I /|2^^2_^2 

A = ^5 ^ P ^. (3.31) 

The quantity A is of order in chiral counting. It is easy to see that, whenever 
the solution of Eq. (j3.3(jp exists, the quantity |Ii| counts at chiral order 1. For 
this reason, the quantity n^^^(i?; q', q) vanishes at 0{p^) in ChPT. 
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4 Equation for the pion-nucleus bound states 

From equations ()2.24|1 and ()2.25j) one can define the pion scattering amplitude 
T{E, ci',q) in the presence of the nucleus 

D{E, q', q) = ^ff^y ~ f + -r, \ ^ TiE, q', q) (4.1) 

which obeys the Lippmann-Schwinger equation 

T{E; q', q) = U{E, q', q) + ^ ^ U{E, q', k) ^^—^^—-^ T [E , k, q) , (4.2) 

with the self-energy Il{E] q', q) playing the role of the potential. Up to and 
including 0{p^), this quantity is given by 

U{E; q', q) = {271 f 5^ {q! - q)n(°)(^' " q') + ^^'\E; q', q) , 

n(i) (E; q', q) = U^'^'^ + H^"") + H^^^) + H^^^) + H^^) , (4.3) 

where the individual terms of Eq. ()4.3|) are given in Eqs. (j3.4p . (j3.17p . ()3.19|) . 
(jT^ and (jSmi), respectively. 

The quantity 11'^°^ (E'^ — q^) defined by Eq. (j3.4|) is both infrared-divergent 
and gauge- dependent. Below we shall demonstrate that, despite this fact, the 
position of the bound-state pole in the scattering amplitude T{E; q', q) is - at 
this order in chiral expansion - both infrared-finite and gauge-independent. In 
order to prove this, we note that, using the theory of scattering on two potentials, 
the scattering amplitude can be given as 

T(£:;q',q) = {2nf5\c^-ciHE'-ci') 

W"{E'-c{ME;^.^W\e^-^'). 

K{s) = [l-A^'^^-{s-ml)B^'\s)]-\ u{s) = {ml-s)[l-K{s)], (4.4) 

where A^^\ B^'^^ are defined in Eq. ()3.4j) . and r(i?; q',q) obeys to the following 
Lippmann-Schwinger equation 

r(E;q',q) = U{E, q^, q) + / U{E, q^ k) r{E, k, q) , 

[/(E,q',k) = K^/^{E^ -c^^)U^^\E;c^,q)K^/\E^ -q^). (4.5) 

From equation ()4.4|) it is clear that the position of the poles in T{E; q', q) and 
t{E] q', q) coincide. On the other hand, the difference between the quantities 
n'^^)(£'; q', q) and [/(£"; q', q) which is due to the factors k^^'^, starts at 0{p^) 
and can be neglected. This is the statement which we aimed to prove: the pole 
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position is determined from the equation ()4.5|) whose potential U{E; q', q) is 
infrared- finite and gauge- independent up to and including 0{p^) in ChPT. Note 
that, since the potential is Hermitian at this order, the position of the bound-state 
pole can also be found by solving the equivalent Klein-Gordon equation 

[E' - q!' - ml]^M) = U{E- q', q) $^(q) , (4.6) 

where $E(q) stands for the wave function of the bound state in the momentum 
representation. 

We finally collect all terms of the potential at this order (replacing cr['^'^l(x) 
by p[°'=^l(x), see Eq. ^^): 

U{E; q', q) = I £^e-^^^'-^>- [U{E; q', q; x) + 0{p')] , (4.7) 

U{E; q', q; x) = - / dh [p%r) + p'{t)] 

J |x — r| 

E - e^F^h 3 , . 4cmlo - 2(C2 + C3)E^ + 2csq! ■ q + Se^F^/i 



2F2 P'^ + ' ^ (^) 



9a_ 
2F2 



E2-2q'-q 3(q'-q)2 -q'2q2i 



+ 



2Mn AMnE^ 



P°(x) 



q'-q , Mr,-Mp. 3 



Equation ()4.8p represents our main result: the complete expression of the optical 
potential for pion scattering on the finite nucleus up to and including 0{p^) in 
ChPT in the presence of electromagnetic interactions and strong isospin-breaking 
effects. Note that, unlike most of previous descriptions, the framework described 
in the present paper enables one to unambiguously obtain the explicit dependence 
of the optical potential on the momenta (off-shell) q and q'. For example, the 
terms proportional to (q' ■ q)^ — q'^q^ vanish in the infinite medium. They 
cannot be obtained from an extrapolation of standard in-medium ChPT or semi- 
phenomenological approaches starting from pion scattering in nuclear matter. 

5 Comparison with existing approaches 

In this section, we compare our results 

i)-iii) to the existing calculations of the in-medium pion mass shift in ChPT j3 

imiia; 

iv) to the pion-nucleus optical potential obtained in Ref. ^2] from ChPT; 
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v) to the empirical pion-nucleus optical potentials which were widely used 
in the literature to describe deeply bound states of pions on the heavy 
nuclei ji]. 

i) In order to compare with the calculations performed in the infinite nuclear 
medium, we put p°(x) = Pp + Pn, P^(x) = Pp - Pn, where Pp,„ = {kp'''^)^ and 
kp'"'^ stands for the Fermi- momentum of the proton (neutron). In the limit of a 
uniform medium, our expression for the self-energy in the absence of electromag- 
netic interactions (e = 0) coincides with the expression for the same quantity at 
0{p^), which is given in Eq. (4.5) of Ref. fT], except for the term proportional to 
g\{Mn — Mp) (or, equivalently, g\c^B{md — rriu) in the chiral expansion). Note 
that, since e = is assumed in Ref. at O(p^) there is no difference whether 
the LEG ci is multiplied by (as in jTT]) or m^o (as in the present paper). For 
the same reason, the contribution involving the LECs /i, /2 is not present in 
Eq. (4.5) of Ref. HI]. 

ii) In Ref. P the threshold self-energy {E = m^) is evaluated at O(p^) in ChPT, 
at vanishing 3-momentum, and with e = 0, rUu = rud, but at kp^ 7^ k^p\ The 
result, given by Eqs. (4)-(5) of this paper, does not depend on the energy E 
as well. Up to and including order p^, this result agrees (apart from the terms 
proportional to and — m„) with our result, if one sets E"^ = m^. However, 
ref. jni does include very important double scattering contributions which first 
appear at 0{jP). 

iii) The results obtained in Ref. ^2] are similar to those of Ref. [H], with the 
exception that in Ref. ^2] the restriction E"^ = is removed. As was mentioned 
above, both jHl determine the pole position from the two-point function of 
the interpolating pion fields. In this quantity one encounters the so-called off- 
shell ambiguity for E'^ 7^ m^, which should disappear in the calculation of the 
pole position in order to be compatible with general principles of quantum field 
theory. In Ref. ^2] it is demonstrated that this actually happens if one expands 
the pole position in powers of the quark mass and eliminates the contributions 
at 0(p'^) and higher (see also [32] )• 

Concerning numerical studies of the in-medium pion mass shift, the perturba- 
tive estimates in all three papers CH [12] yield results around 10 MeV, which 
is too small as compared to the "empirically" deduced shift ~ 25 MeV. In Ref. 
[TT] it is shown, that making the partial resummation of the higher-order terms, 
it is possible to obtain the result for the mass shift which differs from the "per- 
turbative" solution by almost a factor of 2. The main reason for this difference 
is the strong energy-dependence of the self-energy operator. It should be noted, 
of course, that the concept of a pion mass shift in nuclear matter, while of some 
theoretical interest, is of only limited relevance in the context of Coulomb-bound 
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pionic atom states. It is in fact more satisfactory to directly test the theoreti- 
cal predictions for the bound-state observables on the experiment by solving the 
wave equation self-consistently for the finite system - e.g., as done in Ref. [T3] . 



iv) In order to compare with the pion-nucleus optical potential obtained in 
Ref. JH] by extrapolating the results of in-medium ChPT to the finite nuclei, 
we rewrite Eq. fl4.7j) as follows: 

f/(E;q',q;x) = Uc{E;^) + Us{E-^) - {q' ■ q)Up{E-^) 

+ (3(q'-qf -q'V)f/D(i?;x), (5.1) 

where Uc{E]x.) is the Coulomb term (first line of Eq. ()4.8|1 ). and S—,P—,D~ 
wave parts are given by 



2F2 

Acmlo - 2(c2 + C3 - gl/8MN)E^ + 2e^F^f, 



F2 



P°(x) 



p2 r K ; 



E 



p'(x) 



UniE,^) ^ ^^4^P°(x). (5.2) 

The S-wave part of the optical potential coincides with the one of Ref. ^H] 
at 0{p^) in ChPT when setting e = 0. The low-energy constant ci is re- 
lated to the pion-nucleon sigma term, = — 4cim2. The combination 
C2 + C3 receives an important contribution from the A resonance in the 
P-wave ttN amplitude. The values of these LEC's are individually much 
larger that their combination appearing in Us{mT^; x): they must conspire in 
just such a way as to reproduce the observed approximate vanishing of the 
isospin-even nN amplitude at threshold. A small departure from this sub- 
tle balance at threshold is expected to have a large impact on pionic bound 
state energies. This is pointed out in Ref. [12] where it is argued that a 
large part of the "missing S'-wave repulsion" observed in pionic atoms can 
be explained by the strong energy dependence of the low-energy ttN ampli- 
tude in ChPT, while the proper gauge-invariant inclusion of the Coulomb 
potential accounts for most of the remaining effect. This calculation also 
includes the important double scattering terms of 0{p^). The result, con- 
cerning the energy-dependence at 0{p^), is consistent with the findings of 
Ref. pi] for the infinite medium. 

Finally we note that the electromagnetic corrections represented by the 
LECs /i and /2 were so far not taken into account in the literature. 
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• For the P-wave part of the optical potential, the authors of Ref. use 
the time-honored phenomenological parameterization which systematically 
reproduces a large amount of pionic atom data. This parameterization 
includes effects of higher order in density, such as the Lorentz-Lorentz cor- 
rection and pion absorption ^ Ej, which cannot be handled within the 
one-nucleon sector of ChPT. On the other hand, ChPT® makes rigorous 
leading-order statements concerning the P-wave pion self energy for finite 
systems, and it is at this level where a comparison with phenomenology 
is useful. In particular Up{E;x.) can be compared to the leading (linear 
in density) term of the phenomenological parameterization, represented by 
the vrA^ scattering volumes Cq and Ci |2j. (The scattering volumes are 
usually denoted by Cq and Ci, but we need to avoid confusion with the 
LECs Ci, ...,C3.) One may replace E — > in our expression, and rewrite 
the P-wave part of the optical potential as 

^[1 + JJ^jUpi^^-^^) = £oP°(x) - 5,p%^) . (5.3) 

Using the empirical value Cq = 0.21m~^, from Table 6.2 of Ref. |2], fi'yi = 
1.267 and F = 92 A MeV, we obtain cs = -3.2 GeV-^ There is a theoretical 
uncertainty in this number, given by the typical size of next-to-leading order 
corrections. As a consequence, our leading-order estimate of the LEG C3 in 
the low-energy vrA^ effective Lagrangian is not incompatible with the value 
C3 = —4.7 GeV~^ quoted in Ref. [TT]. As for the isovector term, we find 
Ci = 0.17 ai E = m^, in agreement with the value used in the set A 
of Ref. 0. 

Note that the relation between the LEG C3 and the phenomenological P- 
wave scattering volume Cq is different in the case of scattering on a single 
nucleon. Indeed, it is well known that the contribution from the Born 
(nucleon pole) diagrams to this quantity is zero up to and including terms 
of 0{1/M]\[) The differences between the two cases originate in recoil 
effects which differ for a single nucleon and for the nucleus as a whole. In 
the case of scattering on a nucleon the recoil is an effect of order O (p) . In 
the case of a nucleus it is suppressed by an additional factor of 1/A. It is 
therefore natural to find such differences already in our 0{p^) calculation. 
These terms are corrections to the low-density theorem relating the pion 
self-energy to the forward ttN scattering amplitude. It is not possible to 
identify such corrections in the conventional in-medium GhPT, since the 
restriction to p = q makes the separation between P-wave and S'-wave 
ambiguous. 

• There is no D-wave part in the optical potential of Ref. • In our ap- 
proach, this contribution is partly due to the finite-size effect. As already 
mentioned, the presence of such terms cannot be reliably established from 
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the extrapolation of the result obtained at uniform density. They can only 
be controlled in the consistent approach described in this paper. 

v) The phenomenological optical potential for the pion-nucleus interaction close 
to threshold is commonly written in coordinate space as |3] 

2m^t/P^'=° = q{r) + Va(r) V , (5.4) 

where g(r) and a{r) are taken to be energy-independent. From Eq. ()5.2j) we ob- 
serve that there are more terms in the result obtained within ChPT: for example, 
we have in addition the D-wave part of the potential which is absent in Eq. ()5.4|) . 
What is important, the counterparts of q{r), a{r) in ChPT depend also explic- 
itly on the energy E. This energy dependence stems from the underlying chiral 
pion-nucleon dynamics, and is uniquely determined in ChPT. 

A few final remarks are in order. First, we wish to mention that in Refs. jH 
El ^1 ^1 the terms with e 7^ 0, — 7^ in the non-Coulomb part of 
the optical potential have always been neglected. In this paper, we have included 
them on the same footing as all other (strong) terms that arise at the same 
chiral order. Note that the additional electromagnetic terms have sometimes 
led to significant corrections in the observables bound states of hadrons (see 
e.g. miES]). 

As already mentioned, the optical potentials given in Refs. jH [T2] include 
additional terms, non-linear in the baryon density. It is well known that these 
terms give important contributions to the binding energies and widths of pionic 
atoms. In our counting all these terms start at 0{p^) and have not (yet) been 
considered in the present paper. Clearly, in order to achieve a good description 
of pionic atoms, one has to carry out the calculations, outlined in this paper, at 
least up to and including O(p^). 

The focus in this paper is on the presentation of the ChPT® framework for 
finite systems. Numerical algorithms for solving the bound-state equation and 
finding the systematic ChPT expansion of the pole positions are a different issue 
which deserves a separate study. Such an investigation is highly non-trivial. It 
requires checking the convergence of the chiral expansion for the eigenvalues and 
examining the stability of these eigenvalues with respect to small variations of 
the self-energy operator, bearing in mind the strong energy dependence discussed 
previously. 

6 Conclusions 

i) In this paper we propose a novel approach to construct ChPT in the back- 
ground of a finite nucleus: "Chiral Perturbation Theory for Heavy Nuclei" 
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(ChPT® ). It develops the rules for systematically evaluating the pion- 
nucleus optical potential directly for a finite size system. In the present 
paper the full result of the calculations at 0{p^) is presented. 

) The approach is based on certain approximations concerning the description 
of the nucleus containing A nucleons, where A is a large number. The 
content of these approximations, given by Eq. ()2.2H) . is as follows: the 
nucleus in ChPT® does not contain free pions and photons. Secondly, the 
nucleus is not excited by the external pion bound in the atom. 

) For simplicity, it is also assumed that the nucleus is static and in a spin- 
saturated ground state. These assumptions are purely technical and - if 
needed - can be released at a later stage. 

) In ChPT® , the scattering amplitude of the pion in the nuclear background 
is given by a sum of terms, each of which is a product of two factors: the 
scattering amplitude of the pion on 0, 1, 2, ■ ■ ■ nucleons which is systemati- 
cally evaluated in ChPT, and the matrix elements of the normal products of 
the free fermionic fields between the nuclear states. These (and only these) 
matrix elements summarize the necessary nuclear structure information. 

) A crucial ingredient of our approach is the chiral counting for the nuclear 
matrix elements. The bookkeeping based on chiral symmetry allows one 
to carry out a Taylor expansion of these matrix elements. This expansion 
drastically reduces the sensitivity of the final results on the nuclear input. 

) The approach can be pursued systematically to higher chiral orders: 

First, the chiral counting allows one to unambiguously organize different 
contributions, and off-shell ambiguities are absent from the beginning. 

There are no ultraviolet divergences arising at any step, which would signal 
an internal inconsistency: the coefficient functions S**-"^ are finite by con- 
struction, and the phenomenological input on the nuclear matrix elements 
does not introduce ultraviolet divergences either. 

The only general property of the theory which might be lost in the approxi- 
mations that have led to ChPT® , is the gauge invariance of the bound-state 
energy at high chiral orders. We have checked gauge invariance explicitly at 
0{p^). However, at the present stage, we are not aware of a general proof 
that this can be done in all orders. Explicit gauge dependence appearing 
at 0{p^) would signal that, starting from this order, one cannot neglect 
components of the nuclear wave function that contain pions and photons. 
From the phenomenological point of view, the situation is harmless. 
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vii) Our approach is the generahzation of in-medium ChPT which is based on 
the Fermi-gas model. In the hmit of uniform density, our approach reduces 
to the conventional in-medium ChPT. 

viii) In this paper we present a systematic derivation of the leading-order pion- 
nucleus optical potential in ChPT, including isospin-breaking contributions, 
and provide the detailed comparison to other results existing in the litera- 
ture. It is demonstrated that the present approach generates more terms, 
with higher spatial derivatives, than conventional models using the local 
density approximation. In particular, this concerns an additional contribu- 
tion in the P-wave scattering volume, as well as a D-wave term which is 
absent in the conventional parameterizations and can be interpreted as a 
genuine finite-size effect. The negligence of isospin-breaking contributions 
cannot be justified since these come at the same chiral order as the strong 
contributions, and since one anyway includes some of the electromagnetic 
contributions - e.g. the Coulomb potential. An accurate fit to the measured 
bound-state energies of pionic atoms might help to set bounds on the value 
of the electromagnetic low-energy constant /i which enters the expression 
of the optical potential at O(p^). 

ix) Carrying out complete calculations at O(p^) is challenging from several 
points of view. First of all, the contributions at this order, including the 
contributions which are quadratic in the baryon density, are known to be 
important. We do not give any numerical estimates in this paper, based on 
the 0{p^) calculations alone. We plan to address 0{p^) calculations, as well 
as the thorough numerical analysis of the bound-state problem, in future 
publications. 
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Andreas Wirzba for discussions. 

A Chiral Lagrangians 

For convenience, in this appendix we collect chiral Lagrangians, which are used 
in the calculation of the pion self-energy at 0{p^). Notations and conventions are 
identical to the used in Ref. jT2|, from which the formulae below are taken 

+ Cf^ + Cf^ = -\ F,^F^^ - 1 {d,An' 
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7 10 
i=l i=l 

cf^ = E Q^of V , 4'^ = E /.^of . (Al) 

i=l i=l 

Here, ^ denotes the gauge parameter (^ = 1 in the Feynman gauge), and ((■ ■ ■)) 
stands for the trace over the isospin indices. The building blocks for constructing 
the Lagrangian are 

d^U = d^U - iTZ^U + lUC^ , (^^^ ) = i;^ ± + QA^ , 

T^fj.v dfjl^v dyT^ix ^\^fj,i T^u\ ) (^fi^v ^v^ji, ^[-£^5 1 

X = 2B{s + ip) , df,x = df,x - iT^t^X + ix^t^ , ^/.i. = df,A^ - d^A^, , 
C^ = -^[7^„Q], Ci = -i[C„Q], Q = ediag(| -^) , (A2) 
in the meson sector, and 

T/^ = 2 [^^' ^^^] ~ 2 ^^^''^ ~ 2 ^'^'''"^ ' 

x± = ± , x+ = x+-^(x+), 

Q± = ^ (nQnt ± ntQw) , Q± = Q± - ^ (Q±) , Q = e diag(l, 0), (A3) 

in the nucleon sector. Further, i?^, L^, L^,y are defined just like their pionic 
counterparts IZ^, C^, TZ^u, ^^.^i, respectively, with Q replaced by Q. As usual, 
j = {s,p,Vfj_,afj_} denote external scalar, pseudoscalar, vector and axial fields, in 
the form of 2 x 2 matrices. In all expressions we drop the terms which do not 
contain pion fields, and terms of order e^. 

In the above formulae, the symbol e stands for the electric charge. The quanti- 
ties F, m and qa are the pion decay constant, nucleon mass and the nucleon axial 
constant in the chiral limit. The quantity B is related to the quark condensate in 
a standard manner, and the quantity Z is expressed through the charged-neutral 
pion mass difference in the chiral limit 

ml - mlo = 2e'^F^Z + ■■■ (A4) 
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Table 1: Operator basis and the divergent parts (in the Feynman gauge) of the 
LECs in the O(p^) meson Lagrangian [21] and O(e^p^) meson Lagrangian |37j . 
The terms that do not contain pion fields, and terms of order are not displayed. 
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Table 2: Operator basis in the 0{p^) pion-nucleon Lagrangian [2H] and in the 
O(e^) pion-nucleon Lagrangian PSI . 
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Further, the LECs /j and ki are ultraviolet-divergent 
where 

The components of the external sources are defined as 

1 r" r'^ 

^.= 2^J + Y<' «M = Y«^ . = sVr-.", p = rV. (A7) 

In tables ^ and HJ we collect the operator basis for meson and meson-nucleon 
Lagrangians. In addition, in tabled the divergent parts of the LECs Z,, ki are 
listed (see Eq. ^Bj)- 
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